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Abstract. We present Spitzer Space Telescope infrared photometry of a secondary eclipse of the hot Neptune 
GJ436b. The observations were obtained using the 8-/im band of the InfraRed Array Camera (IRAC). The 
data spanning the predicted time of secondary eclipse show a clear flux decrement with the expected shape and 
duration. The observed eclipse depth of 0.58 mmag allows us to estimate a blackbody brightness temperature 
of Tp = 717±35 K at 8 /im . We compare this infrared flux measurement to a model of the planetary thermal 
emission, and show that this model reproduces properly the observed flux decrement. The timing of the secondary 
eclipse confirms the non-zero orbital eccentricity of the planet, while also increasing its precision (e = 0.14 ± 0.01). 
Additional new spectroscopic and photometric observations allow us to estimate the rotational period of the star 
and to assess the potential presence of another planet. 
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1. Introduction 

G J 436 b is one of the few known Neptune-mass extra- 
solar plane ts. It was discover ed by radial-velocity mea- 
surements ([Butler et al.l 120041) as a planet with a period 
of 2.6 days and a minimum mass of 21 M®. Follow-up 
Doppler observations of GJ 436 refined the planetary mass 
and the or bital parameters, in cluding an eccentricity of 
0.16 ±0.02 (iManess et al.ll2007l hereafter MOT). Our team 
( Gillon et al.ll2007al hereafter G07a) discovered the tran- 
siting nature of GJ436b, enabling us to measure a plane- 
tary radius 4 . This discovery and the corresponding 
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measurements of the planetary radius and mass indicated 
a planet composed mostly of ice, probably surrounded by 
a small H/He envelope. 

Because of the small size of the parent star {R ~ 
0.4 Rq) and the short orbital period of GJ436b, the 
planet-to-star luminosity ratio in the infrared is compa- 
rable to that of many known hot Jupiters, despite the 
planet's much smaller radius. Furthermore, the M dwarf 
GJ436 is rather bright in the infrared (K ~ 6). Detection 
of the thermal emission from this small planet had thus 
been expected to be within the reach of the Spitzer Space 
Telescope. 

Following our transit discovery, we submitted a 
Discretionary Director Time (DDT) Spitzer proposal to 
better characterize this interesting planet. We applied 
for photometric observations of the primary transit us- 
ing the 8-/xm band of the InfraRed Array Camera IRAC 
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(| Fazio et al.l 12004 ) in order to get a very accurate ra- 
dius measurement and constrain the bulk composition 
of the planet. We also applied for photometric obser- 
vations of the secondary eclipse in the four bands of 
IRAC (3.6, 4.5, 5.8 and 8 /x m), in the 16-um ba nd of the 
InfraRed Spectrograph IRS (jHouck et al.l |2004[) and the 
24- /xm band of the Multiband Imaging Photometer MIPS 
(jRieke et al.l 120041 ) to assess the atmospheric tempera- 
ture, albedo, heat distribution efficiency and composition. 
However, the observations were actually triggered and per- 
formed as part of an existing Target of Opportunity (ToO) 
program (ID 30129, PI J. Harrington) which has a to- 
tal priority for the observations of transiting planets. The 
main goal of this ToO is to deliver to the community with- 
out any proprietary period optimal Spitzer observations 
of transiting planets. 

Spitzer observed the transit and the secondary eclipse 
of GJ436 in the 8-/xm IRAC band on June 29 and 30 
respectively. The data of the primary transit were made 
publicly available on July 13th 2007. Our direct analysis of 
these data allowed us to det ermine a very accur ate radius 
for G J 436b (Rp = 4.2 Rm. lGillon et"ani2007bl . hereafter 



GOTb) and to confirm the presence of an H/He envelope. 
Spitzer data of the secondary eclipse were not released to 
the community until July 17th 2007, due to an oversight 
that occurred at the Spitzer Science Center. This explains 
why we separated our analysis and present here our results 
regarding the secondary eclipse data. 

During the wri ting of this study, a paper by 
Deming et al.l ( 2007 ) reporting primary and secondary 
eclipses analyses has been submitted to ApJ and put on 
astro-ph. The present analysis has been conducted inde- 
pendently from their work. Their results are consistent 
with the ones presented here. 

Analyzing the secondary eclipse data, we report here 
the detection of a secondary eclipse and draw conclu- 
sions about the thermal emission of G J 436 b and refine 
its orbital parameters, allowing a better understanding of 
GJ 436 dynamics by exploring the contingency of a sup- 
plementary planet. 

In addition, we report here on additional ground based 
observations to determine the stellar rotational period. We 
followed the photometric intensity and the Ca II H-l-K 
activity index of GJ 436. Although the photometric data 
are sparse and cover only 50 days, we find some evidence 
that the stellar rotational period is of the order of 50 
days, which is also consistent with long-term Call mea- 
surements. 

Section 2 describes the observations and the reduc- 
tion procedure. Our analysis of the obtained secondary 
eclipse time series is described in Section 3. In Section 4, 
we analyze the infrared emission from the planet and draw 
some conclusions about its atmosphere composition. We 
detail an orbital analysis, encompassing the possibility of 
a perturbing planet, stellar activity and GJ436b orbital 
parameters refinements in Section 5. Our conclusions are 
presented in Section 6. 



2. Observations and data reduction 

2.1. Spitzer IRAC observations 

GJ436 has been observed on June 30th UT for 6 hours, 
to cover the secondary eclipse, resulting in 49920 frames. 
Observations were made so as to encompass the expected 
secondary eclipse window, whose timing calculations were 
made by taking into account transit timing and orbital 
eccentricity. Due to the uncertainties on eccentricity and 
argument of periastron, a larger time- window was chosen 
to ensure the detection of the secondary eclipse. Data ac- 
quisition was made using IRAC in its 8-/im band with the 
same mode and strategy employed for the primary transit 
(G07b). 

We combine each set of 64 images using a 3-cr clipping 
to get rid off transient events in the pixel grid, yielding 780 
stacked images for the secondary eclipse, with a temporal 
sampling of ~ 28s. Heliocentric Julian Day (HJD) con- 
version was made according to the mean Spitzer orbital 
position at the time of each exposure and GJ436 appar- 
ent position. Spitzer position ephem erides were obtained 
through JPL-Horizons web interface ( Giorgini et al.l[r996h 
and converted from TT (Terrestrial Dynamic Time) to 
UTC. 

We faced the same instrumental rise issue noticed in 
our work on primary transit. To mitigate its effect, we zero 
weight the eclipse and the first 100 points of the time- 
series. We then divide the lightcurve by the best fitting 
asymptotic function with three free parameters and eval- 
uate the average flux outside the eclipse to normalize the 
time series, exactly as for the primary transit. The rms of 
the resulting time series evaluated outside the eclipse is 
the same as for the primary (G07b): 0.7 mniag, which is 
1.2 times GJ436's photon noise. 



2.2. Ground-based photometry 

To assess the variability of the star, we observed GJ 436 
with the Euler Swiss telescope located at La Silla 
Observatory (Chile) and the Frangois-Xavier Bagnoud 
Observatory's (OFXB) 0.6m telescope located at Saint- 
Luc (Switzerland). Observations occurred in 14 nights 
from May 4th to May 21th. A sequence of 10 exposures 
was done every night. The same strategy used for our 
observation of the May 2nd transit (G07a) was applied 
(V-band filter, 80s exposure time, defocus to ^ 9"). The 
data reduction was also similar. We also use for our anal- 
ysis of the GJ 436 variability the May 2nd out-of-transit 
data and the photometric lightcurves obtained with the 
OFXB 0.6m telescope during our search for the transits of 
G J 436 b (G07a). We scale OFXB points with Euler ones 
because of the filters slightly different bandpasses. At the 
end, our data amounts to 24 points spanning 48 days. The 
lightcurve is represented in Fig. 6, and discussed in Sect. 
5.3. 
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2.3. Ground-based spectroscopy 

Since the discovery of G J 436 b (|Butler et al.ll2004l ) , we ob- 
tained additio nal spectra of the s tar with the ESO Harps 



spectrograph (jMavor et al.l 120031 ). Harps is mounted on 



ESO 3.6m telescope and is dedicated to high precision 
radial-velocity measurements thanks to its resolution of 
IIO'OOO and a wavelength range coverage between 3800 
and 68OOA. To assess the stellar activity and rotation we 
used 23 high SNR spectra from which we measured the 
Ca n H-l-K index. Results are discussed in Sect. 5.3 

3. Analysis of secondary eclipse time series 

We fit a non-limb-dark ened eclipse profi l e to th e secondary 
eclipse data using the iMandel fc Agoll (j2002r ) algorithm. 
The eccentricity of the orbit is considered as described 
in G07b, taking the values for the eccentricity e and the 
argument of periastron uj from M07. The formula connect- 
ing w to the true anomaly / at the orbital location of the 
secondary eclipse is: 



/ 



2 



- ui. 



(1) 



We fix the stellar and orbital parameters to the values 
mentioned in G07a. The free parameters are the central 
epoch of the secondary eclipse Ts and the flux decrement 
AFs- The fit procedure and the error bars estimation is 
similar to the one described in G07b. The obtained value 
for Ts and AFg, including their respective error bars are 
given in Table 1. Figure 1 shows the best-fit theoretical 
curve superimposed on the lightcurve (zoomed on sec- 
ondary eclipse center, binned for clarity) and the residuals 
of the fit. 

After having derived an accurate value for the eccen- 
tricity (see Section 5), we perform a new fit to the sec- 
ondary eclipse, taking into account the new values for the 
orbital eccentricity and the true anomaly at the orbital 
location of the eclipse, and their new error bars. The ob- 
tained values are in excellent agreement with the one given 
in Table 1. 

4. Infrared radiation 

While GJ436b is properly classified as a hot Neptune, 
the irradiation by the host star is weaker than for most 
hot Jupiters. Consequently, the contrast measurement re- 
ported here is that of the coolest exoplanet atmosphere 
detected so far. The atmospheric temperatures are pre- 
dicted to be low enough for carbon to be bound in CH4 
(instead of CO as is the case for most hot Jupiters), plac- 
ing G J 436 b in a yet unexplored exoplanet atmospheric 
regime. The temperatures should also be cool enough for 
NH3 absorption to appear between 10 and 11 /im. This 
situation is comparable to T dwarfs, which have promi- 
nent absorption bands of NH3 at 10.5 /xm as seen in re 
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Fig. 1. Top: Zoomed binned time series for the secondary 
eclipse. The best-fit theoretical curve is superimposed. 
Although unbinned data were used for the fit, points are 
binned by 5 for plotting purposes. Bottom: The unbinned 
residuals of the fit. Their rms is 0.7 mmag. 
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Fig. 2. Model planet-star flux ratios for GJ436b assum- 
ing that the absorbed stellar flux is redistributed across 
the dayside only (top curve) and uniformly redistributed 
across the entire planetary atmosphere (lower curve). In 
both models the composition is equal to that of the host 
star. For the wavelength range shown, the majority of the 
planet spectral features are produced by water, methane, 
and ammonia absorption. The filled black diamond is our 
Spitzer contrast measurement at 8 /im with its associ- 
ated error bars while white diamonds are the model con- 
trast values in the Spitzer IRAC and IRS bandpasses. 
The dashed line is the contrast curve for 700K blackbody 
planet spectrum. 



cent Spitzer IRS observations (jCushing et al 
Fig. 2 we compare our 8 /um contrast measurement to 



synthetic planet-star flux ratios calculated follow ing the 
methods described in Barman et al.l ( 2001 . I2OO5I) for two 
different assumptions for the day-to-night energy redis- 
tribution. The hotter dayside model corresponds to no 
20061 ). In redistribution of energy to the night side, while the sec- 
ond (lower flux) model assumes very efficient redistribu- 
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Mid-SE timing [HJD] 
Flux decrement [AFs] 
Tb at 8 /im [K] 
Orbital eccentricity 



2454282.333 ±0.001 
0.00054 ±0.00007 
717 ±35 
0.14±0.01 



Table 1. Parameters derived from the secondary eclipse 
for GJ436b. SE stands for Secondary Eclipse. 



tion of energy capable of completely homogenizing the day 
and night sides. As can be seen, our 8 /im measurement 
agrees very well with the hotter of the two models sug- 
gesting that redistribution is fairly inefficient. However, 
it is impossible to constrain the bolometric flux emerg- 
ing from the planet (and thus the true energy budget of 
the day and night sides) with a single flux measurement 
in one bandpass. If energy redistribution is highly depth- 
dependent, as indicated by recent dynamical simulations 
( Cooper fc ShowmanI 2005[ l. then it remains possible that 
significant amounts of energy is being transported to the 
nightside, resulting in a warm nightside and cooler day- 
side at depths above or below the 8 /im photosphere. The 
agreement with the model spectrum suggests that obser- 
vations at other Spitzer bandpasses should be possible and 
will allow further valuable constraints on both the atmo- 
spheric composition and the energy redistribution. In par- 
ticular, the TOOK blackbody planet spectrum (dashed line, 
Fig. 2) illustrates the value of observations at 4.5 and 16 
/im as helpful probes of different atmospherics depths hav- 
ing different brightness temperatures. Here, we estimate 
a temperature of Tb — 717 ±35 K at 8 /im, by compar- 
ing the observed contrast to blackbody SEDs divided by a 
synthetic stellar spectrum (Teff — 3350 K, M07), weighted 
by the radii ratio squared. We then varied the blackbody 
temperature until the 8 /im integrated contrast matched 
the observed contrast value. 



5. Orbital analysis 

5.1. The non-zero eccentricity 

One noticeable characteristic of G J 436 b is its non-zero 
eccentricity (e=0.16±0.02 - M07). It contrasts with most 
known short-period exoplanets (P < 5 days) which have 
very small eccentricities, often indistinguishable from zero. 
Unfortunately, moderate eccentricities are difficult to con- 
strain with radial- velocity measurements, and M07 warn 
that the quoted errors of the orbital parameters, based 
on the bootstrap technique, may lead to wrong estimates 
in some cases. To assess the statistical significance, they 
choose to use a rigorous Bayesian analysis and found the 
eccentricity to be greater than with a high confidence 
level. Still, GJ436b's eccentricity is only known with a 
large uncertainty. 

To improve the determination of GJ 436 b's eccentric- 
ity we combine Spitzer eclipse timings with M07 radial 
velocities and perform a combined fit. As M07 have shown 
with a high confidence level, a positive radial-velocity 
trend is present in their data, we choose a model made 



of a planet plus a linear drift. Our minimization is based 
on the Levenberg-Marquardt algorithm (Press et al. 1992) 
and, as a maximum likelihood approximation, minimizes 
the following x^: 



Tp — Tp 



Ts-Ts 



(2) 



where Vi is the i radial velocity given in M07 and Tp 
and Ts are respectively the timings of the Spitzer primary 
transit and secondary eclipse reported in this paper. The 
corresponding error estimates are e„,i, erp and cts, and 
Vi, Tp and Ts are their corresponding computed value, 
according to the chosen model. 

We find the to be minimum with an orbital period 
P = 2.643859 days, a semi-amplitude K = 18.2 ms~^, a 
date of the passage at periastron Tq — 2454198.2056714 
HJD, an argument of periastron lu — 350°, an orbital 
eccentricity e = 0.14±0.01, a radial- velocity offset 7 = 
4.2 ms~^ and slope dv/dt = 1.4 ms~"'^yr~"'^. 

For this fit, the squared root of the reduced is 1.84, 
marginally higher than a fit with radial velocities alone 

( V X \rv only ~ 1-81). 

To derive the error of fitted orbital parameters we sim- 
ulate 1000 virtual sets of new radial velocities and new 
eclipses timings. In each set, the radial -velocity data are 
randomized with a bootstrap algorithm ( Press et al.lll99sj) 
and the eclipses timings are randomly generated according 
to a normal distribution, with mean and standard devia- 
tion given by the actual timing values and their error, re- 
spectively. Figure 3 shows these probability distributions 
for the eccentricity in both cases, when only the radial- 
velocity data are used and when a combined fit of radial- 
velocity and eclipses timings data is performed. The de- 
termination of eccentricity is clearly improved by the ad- 
dition of eclipses timings, which bring the l-cr error on e 
down to 0.01. 

Spitzer observations therefore strongly confirmed 
G J 436 b unusual eccentricity. M07 pointed out that it may 
be due either to its own structure (i.e. a high tidal-quality 
factor Q) or to an additional long-period companion peri- 
odically interacting with the planet and pumping up its ec- 
centricity. GJ 436 b has since been caught in transit and we 
now have a precise measurement of its radius. Considering 
GJ 436 is probably more than few billion years old, we can 
estimate what Q would dissipate the tidal circularization 
up to this age. 

T o match an age >2 Gy r, a. Q > 10^ is neces- 
sary ( Adams fc Laughlinll2006r ). which is much more than 



Neptu ne in the solar system for which iBanfield fc Murray 



19921) give 1.210^ < Q < 3.310^. Thus, interaction 
with another companion is the most likely explanation for 
GJ 436 b's large eccentricity, probably due to the long pe- 
riod companion suspected from the radial-velocity trend 
in M07 data. 
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Fig. 3. Probability distributions for the eccentricity re- 
sulting from randomly generated datasets including: Top: 
Radial velocity data only. Bottom: Radial velocities + 
transit and secondary eclipse timings. 



5.2. Looking for additional planets 

The improvement in the determination of orbital parame- 
ters provides an opportunity to look for additional planets 
in the radial- velocity data. Such analysis is also motivated 
by the ^/X^ of our solution, which is larger than one. 

A period analysis of the residuals around the best so- 
lution (Fig. 4) shows no significant power excess at any 
period. The highest peak is found at P '--^ 5.602 days and 
is attributed a 92% false alarm probability by bootstrap 
randomizations. In conclusion, except the companion sus- 
pected in Sect. 5.1, the present data set shows no evidence 
for additional low mass exoplanets in the GJ 436 system. 

5.3. Investigating residuals: the stellar activity 

An alternative way to explain that the dispersion of the 
radial-velocity residuals is in excess compared to the in- 
ternal errors is to invoke the stellar activity. If present 
on the stellar surface, spots are known to modulate the 
Doppler measurement and to introduce 'jitter' or ad- 
ditional coherent signa l in radial-velocity measurement 
( Saar fc Donahudll997t) . 

Earlier this year we published the discovery of a 
m sin i = 11 Mm planet o rbiting the nearby M dwarf 




51500 52000 52500 53000 53500 54000 
Juhan date -2,400,000 [day] 



.03. 
P 

cn 



40 



20 




10 100 
Period [day] 



1000 



Fig. 4. Residuals around the best-fit orbital solution. Top: 
0-C with their error bars. Bottom: periodogram computed 
from residuals. 
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GJ674 (jBonfils et al.ll2007^ . In addition to the Doppler 
signal induced by the planet, we clearly identified a sec- 



Fig. 5. Ca II H+K periodogram obtained from high SNR 
spectra with Harps spectrograph. The arrow points at 
the power excess around P ~ 48 days. 



ond signal of period ~ 35 days in the residuals of the 
one-planet fit. We have shown that Ca II H&K emission 
lines were varying in phase with this second signal, demon- 
strating it was due to a spot rather than a planet. This 
analysis was given further credit by a clear photometric 
counterpart to the spectral-index variation. 

To investigate the activity of GJ 436, we can thus apply 
the same spectroscopic diagnostic as we did for GJ674, 
thanks to Harps spectra we obtained since 2004. Figure 
5 hence represents the periodogram of Ca II H+K index 
measured on 23 high SNR spectra of GJ436. It displays 
a power excess around P ^ 48 days that identifies the 
rotation period of GJ 436. Bootstrap randomizations give 
a false alarm probability < 1% for this peak. 

Moreover, complementary photometric observations 
we did to monitor the long-term activity of GJ436 (Fig. 
6) confirm that a spot is present on GJ 436 surface and 
that the rotational period is likely more than 40 days. On 
a 50 day-time span the variation of the flux has an ampli- 
tude of ~ 1%. We know from the spectral index variation 
that 50 days is close to the rotational period and it is thus 
reasonable to assume this amplitude for the photometric 
signal. With an estimate of the amplitude of the photo- 
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Fig. 6. Long-term lightcurve obtained with the Eulcr 1.2m 
telescope at ESO La SiUa Observatory and the 0.6m tele- 
scope at FXB Observatory. 

metric variation, plus an approximate rotational period, it 
becomes possible to estimate the amplitude of the activity 
in duced radial- veloc it y vari ation. 

baar fc Donahue! (|l997h have done some simulations 



and found the radial-velocity amplitude Kg induced by a 
spot follow approximately the relation : 



Ks\t[1S 



6.5 X /"-^ X usini. 



(3) 



where fs is the size of the spot (expressed in percent of the 
stellar disk) and usini is the projected rotational velocity 
of the star. 

In the case of GJ674, considering its radius (0.34 Rq) 
and its rotational period (34.8 days), we calculate a v sini 
of 0.5 kms^^. Equation 3 then converts the observed flux 
variation (~ 2.6%) into a radial-velocity amplitude Kg ^ 
8 ms~^, close to the measured amplitude (6 ms~^). 

The same numerical application for GJ436, with a ra- 
dius Ri, — 0.463 R© (G07b), a rotational period Prot ^ 45 
days, and a filhng factor fs^l% lead to Kg ^ 3 ms^^. 
The spot is thus responsible for a typical dispersion of 
~ 2 ms^^, which, co-added to the typical radial- velocity 
errors (~ 2.4 ms^^), explains most (if not all) of the dis- 
persion observed for the residuals around our best solu- 
tion (~ 4 ms""'^). Ultimately, to better weight the errors 
between radial- velocity data and eclipses timings data, we 
introduce this 'jitter' in our fitting procedure. Its impact is 
negligible as the estimated parameters remain unchanged. 

6. Conclusions 

Since the discovery, G J 436 b has showed itself as a pecu- 
liar planet and has risen a strong interest from the commu- 
nity regarding its composition or supplementary planets 
in the system. Spitzer data gathered from the primary 
and secondary eclipse are of great help to answer some of 
those questions as discussed in G07b and in this present 
study. 

We especially learn from the infrared emission mea- 
surements at 8 /xm and planetary atmospheres models 



that G J 436 b is characterized by an envelope composed 
of H, He, H2O and CII4. Also, our contrast measurement 
is consistent with a model planet that has very inefficient 
day-to-night redistribution at 8 /im photospheric depths 
on the dayside. 

Moreover, transit and secondary eclipse respective 
timings combined with radial velocities prove that 
G J 436 b has an eccentricity significantly greater than 
zero. Gonsidering a reasonable tidal dissipative factor, we 
estimated the orbital circularization timescale to be likely 
shorter than GJ436 age. We therefore conclude that the 
non-zero eccentricity is probably the result of a dynami- 
cal interaction with an additional companion in the sys- 
tem, maybe the long period companion suspected from 
the radial-velocity trend in MOT data. 

In the course of our orbital analysis we try to find an 
additional planet around GJ436, but no significant peri- 
odicity is found in the residuals of our best fit. Gonversely, 
we identify that GJ 436 has a spotted surface and proba- 
bly rotates with a period Prot ~48 days. We estimate that 
this magnetic activity noises the radial-velocity signal at 
a level of ~ 2 ms~^, therefore explaining most (if not all) 
the residual dispersion around our best solution. 

Nevertheless, the full potential of Spitzer concerning 
G J 436 b has not been explored yet, especially regarding 
thermal emission spectral coverage. Complementary ob- 
servations in the 3.6, 4.5, 5.8-/xm IRAC, 16-/im IRS and 
24- /xm MIPS channels are due between Nov. 2007 and Feb. 
2008. They will certainly bring new constraints on the at- 
mosphere composition of this planet. 
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